Majority- Vote Model on a Random Lattice 
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The stationary critical properties of the isotropic majority vote model on random lattices with 
quenched connectivity disorder are calculated by using Monte Carlo simulations and finite size 
analysis. The critical exponents 7 and f3 are found to be different from those of the Ising and majority 
vote on the square lattice model and the critical noise parameter is found to be q c = 0.1 f7 ± 0.005. 
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I. INTRODUCTION 

Random lattices play an important role in the descrip- 
tion of idealized statistical geometry in a great variety 
of fields 0, 0, 0, Besides its potential to describe 
the topology in several condensed matter problems, the 
comparison of the universality class of systems between 
random and regular lattices is also the subject of in- 
tense research IE IS HI- The randomness consid- 
ered in these studies either by random variation of the 
coupling strengths, random deletion of bonds or sites, 
or by spatially correlated lattices. In particular, accord- 
ing to the Harris criterion |jj valid for the random-bond 
paradigm, random disorder is marginally important to 
the two-dimensional Ising model since the specific heat 
critical exponent is a — 0. However, this criterion could 
not be applied to lattices with a non-periodic coordina- 
tion number, for these models, Luck[To| formulated a 
criterion to such cases. For example, Jankc ct al. flU l us - 
ing the single-cluster Monte Carlo update algorithm O , 
reweighting techniques |l3j ] and size scaling analysis 3| j 
simulate the Ising model in a two-dimensional Voronoy- 
Delaunay random lattice. They calculated the critical 
exponents and found that this random system belongs to 
the same universality class as the pure-two dimensional 
ferromagnetic Ising model. 

The Voronoi-Delaunay network is a type of disordered 
lattice exhibiting a random coordination number that 
varies from 3 to 00, depending on the number of sites. 
In addition, the distance r between nearest neighbors 
changes in a random way from site to site. Lima et 
al. used a two-dimensional Voronoi-Delaunay lattice to 
study the ferromagnetic Ising model and the Potts 
model [l6j. As the bond length between the first neigh- 
bors varies randomly from neighbor to neighbor, they 
considered a coupling factor decaying exponentially as: 




J, 



J exp(- 



(1) 



where ry is the relative distance between sites i and 7, 
Jo is a constant, and a > is a model parameter. In [la ] 
they calculated the critical point exponents 7/V, (3/v and 
v, and concluded that this random system belongs to the 
same universality class as the pure-two dimensional fer- 
romagnetic Ising model. In they studied the three- 
state Potts model and found that critical exponents 7 



FIG. 1: Magnetization m(q) as a function of the noise param- 
eter q for several values of the system size N. 



and v are different from the respective exponents of the 
three-state Potts model on a regular square lattice. How- 
ever, a ratio 7/1/ remains essentially the same. They also 
found numerical evidences that the specific heat on this 
random system behave as a power-law for a — and as 
a logarithmic divergence for a = 0.5 and a — 1.0. 

It has been argued that nonequilibrium stochastic spin 
systems on regular square lattice with up-down symme- 
try fall in the universality class of the equilibrium Ising 
model 0- This conjecture was found in several models 
that do not obey detailed balance 0, 0, |2(J ■ Campos 
et al. |2l| investigated the majority-vote model on small- 
word network by rewiring the two-dimensional square lat- 
tice. These small- world networks, aside from presenting 
quenched disorder, also possess long-range interactions. 
They found that the critical exponents 7/1/ and (3/v are 
different from the Ising model and depend on the rewiring 
probability. However, it was not evident that the ex- 
ponents change was due to the disordered nature of the 
network or due to the presence of long-range interactions. 

Here, we analyze a relatively simple nonequilibrium 
model with up-down symmetry in a random lattice with 
quenched conectivity disord er ( a = 0), namely the 
isotropic majority vote model |22j j. Our main motivation 
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FIG. 2: Log- log plot of the magnetization at q = q c versus L. 
The solid line is the best fit with slope —/3/u = —0.112(4). 

is to investigate whether only the presence of quenched 
lattice disorder is capable of modifyng the exponents 7 jv 
and /3/v using the Voronoi-Delaunay random lattice. Our 
numerical results suggest that the critical exponents, in 
the stationary state, are different from those of the Ising 
model and the isotropic majority vote model on a square 
lattice. In what follow we will utilize Monte Carlo simu- 
lations and finite-size analysis. 

II. MODEL AND SIMULATION 

For each point in a given set of points in a plane, we 
determine the polygonal cell that contains the region of 
space nearest to that point than any other. We consid- 
ered two cells neighbors when they possess an extremity 
in common. From this Voronoi tessellation, we can ob- 
tain the dual lattice by the following procedure. 

(a) when two cells are neighbors, a link is placed be- 
tween the two points located in the cells. 

(b) From the links, we obtain the triangulation of space 
that is called the Delaunay lattice. 

(c) The Delaunay lattice is dual to the Voronoi tessella- 
tion in the sense that its points correspond to cells, links 
to edges and triangles to the vertices of the Voronoi tes- 
sellation. We consider a two-dimensional Majority Vote 
model, on this Poissonian random lattice, defined |2"3,l23j 
by a set of "voters" or spins variables {o^} taking the val- 
ues +1 or —1, situated on every site of a Delaunay lat- 
tice with N sites and periodic boundary conditions, and 
evolving in time by single spin-flip like dynamics with a 
probability Wi given by 

Wi{(l) = -[1 - (1- 2q)a l S(J2^+ delta)} (2) 

where S(x) — sign(x) if x ^ 0, S(x) = if x = 0, and 
the sum runs over all nearest neighbors of <Tj. In this 
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FIG. 3: Reduced fourth-order cumulant U(q) as a function of 
q for several values of JV. Within the accuracy of our data, 
all curves intersect at q c = 0.117(5). The value of U(q) at the 
intersection is U* = 0.61(2). 

lattice, the coordination number varies locally between 3 
and co. The control parameter q plays the hole of tem- 
perature in equilibrium systems and measures the prob- 
ability of aligning antiparallel to the majority of neigh- 
bors. For simplicity, the length of the system is defined 
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FIG. 4: Log-log plot of dU(q)/dq at q — q c versus L. The 
solid line is the best fit with slope \ jv — 0.94(6). 

here in terms of the size of a regular lattice, L — N 1 / 2 . 
We perform simulations for different lattice sizes JV = 2 4 , 
where i varies from 8 to 14. For each size, we gener- 
ated 50 randomly chosen lattice realization, where each 
simulation started with a random configuration of spins. 
From a given configuration, the next one was obtained 
as follows: (a) Choose a spin {i} at random, (b) Gen- 
erate a random number r uniformly distributed between 
zero and one. (c) Flip spin i when r < W{(a). In our 
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simulations, 5xl0 4 Monte Carlo steps were required for 
attaining the stationary state. After that, we estimated 
the averages, for any lattice size, using 10 5 Monte Carlo 
steps. 

We define the variable m = Cj/iV. In particular, 

we were interested in the magnetization [23, |24( , suscep- 
tibility and the reduced fourth-order cumulant: 



m (q) = [< \ m \ >}av 



X(q) = N[< m 2 > - < |m| >%, 



(3) 



(4) 
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where < ... > stands for a thermodynamics average and 
[...]„„ square brackets for averages over the 50 realiza- 
tions. We calculated the error bars from the fluctuations 
among realizations. Note that these errors contain both, 
the average thermodynamic error for a given realization 
and the theoretical variance for infinitely accurate ther- 
modynamic averages which is caused by the variation of 
the quenched random geometry of the 50 lattices. 

These quantities are functions of the noise parameter 
q and obey the finite-size scaling relations 

[< \m\ >} av = L-^f m (x)[l + ...}, (6) 
X = L-^f x (x)[l + ...}, (7) 

^-=L^f u (x)[l + ...}, (8) 
dq 

where v, [3, and 7, are the usual critical exponents, fi(x) 
are the finite size scaling functions with 

X = (q - q c )lM v (9) 

being the scaling variable, and the brackets [1 + ...] indi- 
cate corrections-to-scaling terms. 

III. RESULTS 

In Fig. 1 the magnetization is shown as a function of 
the noise parameter (q) for several values of L. As can be 
noticed there is a phase transition from an ordered state 
(Ml > 0) to a disordered state (Ml ~ 0). This figure 
displays that for q > q c the magnetization disappears 
when larger values of L are considered, whereas it reaches 
a finite value for q < q c . In Fig. 2 a log-log plot of the 
magnetization at q — q c versus L is shown. The ratio 
between the critical exponent (3/v = 0.112 ± 0.004 is the 
slope of the straight line fitted to the data points. Within 
the numerical accuracy, we found that these exponents 
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FIG. 5: Susceptibility xil) as function of q for several values 
of N. 



are distinct from the exponents characterizin g th e class 
of universality of the equilibrium Ising model 22] . 

To determinate the critical point the reduced fourth- 
order cumulant U is plotted against q. In Fig. 3, the 
critical point q c is estimated when all curves, for dif- 
ferent sizes L, intersect in the same point. We get 
q c = 0.117 ± 0.005 and U* = 0.61 ± 0.02. The value 
of U* is, considering the error bar, the same as the one 
obtained for Ising model on a square lattice with periodic 
boundary conditions U* = 0.611 ±0.001. Another way of 
getting the critical point is through the relation of scale 
q^ax( L ) ^qc + aL- 1 ^, where we get q c = 0.117 ±0.003 
using v — 1.06. 

To obtain the critical exponent ^, we calculated nu- 
merically U (q) = dll(q)/dq at the critical point for each 
value of L. The results are in good agreement with the 
scaling relation (7). Then, we plotted lull versus InL, as 
displayed in Fig 4. The straight line represents the best 
fit to the data points. The slope gives 1/v = 0.94 ±0.06, 
which corresponds to v — 1.06 ± 0.08. 

In Fig. 5, we have ln\(q) as a function of q for several 
values of L. In order to study the universality of the 
model, the ratio 7/zv was estimated from the log- log plot 
of the value of the susceptibility x(lc) versus L. We 
estimated the critical exponent 7/1/ = 1.51 ± 0.04 from 
the best fit of data points, as displayed in Fig. 6. 



IV. CONCLUSION 

We have studied the majority- vote model on Voronoi- 
Delaunay random lattices with periodic boundary con- 
ditions. These lattices possess natural quenched dis- 
order in their connections. We verified whether only 
this type of disorder is significant to obtain critical ex- 
ponents different of those found for this model in the 
regular lattice (that have the same exponents of the 
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FIG. 6: Log-log plot of the susceptibility at q c versus L. From 
it we estimate the critical exponent y/u = 1.51(4) as the best 
fit of data points. 

Ising model in two dimensions). We measure the ex- 
ponents 7/V, j3Jv and v. The best fit of these expo- 



nents provided v = 1.06 ± 0.08, 7/1/ = 1.51 ± 0.04, and 
f3/i> = 0.112 ± 0.004 and U* = 0.61 ± 0.02. The criti- 
cal exponents (3/v and "f/v are different from the exact 
values of the Ising model and majority-vote model on a 
regular square lattice p2l |. Our results are in a gree ment 
with the results obtained by Campos et al. [21j that 
studied this same model on a small world network pre- 
senting quenched disorder and long-range interactions. 
They found that the critical exponents depends on the 
shortcuts introduced in the network. In summary we 
showed here that the presence of quenched connectivity 
disorder is enough to alter the exponents (3/v and 7/1/ of 
the pure model and therefore that disorder is a relevant 
term to such non-equilibrium phase-transition. However, 
the critical value of the forth-order cumulant remains the 
same as that of the pure model. 
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